Is Sr 2 Ru04 a triplet superconductor? 
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The field dependence of the specific heat "/(H) at lower temperatures in Sr2RuC>4 is analyzed 
by solving microscopic Eilenberger equation numerically. We find that systematic y(H) behaviors 
from a concaved yH to a convex H a (a > 1) under H orientation change are understood by taking 
account of the Pauli paramagnetic effect. The magnetizations are shown to be consistent with it. 
This implies either a singlet pairing or a triplet one with d-vector locked in the basal plane, which 
allows us to explain other mysteries of this compound in a consistent way. 

PACS numbers: 74.20.Rp, 74.70.Pq, 74.25.Op, 74.25.Bt 



Superconductors are classified into two distinctive 
groups, either spin-singlet or spin-triple pairings. While 
almost all superconductors, including high T c cuprates 
belong to the former, the later is extremely rare and dif- 
ficult to find. Only a few examples of superconductors 
are discussed for its possibility; In a heavy Fermion ma- 
terial UPt3 the identification of a triplet pairing has been 
firmly established P, Q • The observed multiple phase di- 
agram in field (H) versus temperature (T) plane, consist- 
ing of three phases A, B and C, is reasonably explained 
only in terms of triplet pairing. This situation is sim- 
ilar to superfluid 3 He where two subphases ABM and 
BW are identified in pressure vs. T plane 0]. Knight 
shift (KS) experiment by NMR has played a fundamen- 
tal role to confirm the theoretical predictions in UPt3. It 
was particularly crucial that both field directions where 
KS is changed and unchanged below T c are found ex- 
perimentally [1] as predicted 0], identifying the d- vector 
direction. 

Sr2Ru04 is second prime candidate for a triplet pairing 
superconductor [5j . A variety of theoretical and exper- 
imental works have been devoted to establishing it, but 
it turns out after a decade of its discovery [6| that it 
is extremely difficult to identify the spin structure of a 
Cooper pair although the gap structure with line node is 
well established now. For example, it is pointed out that 
recent phase-sensitive experiments by Nelson et al. 0], 
Kidwingira et al. Q and Xia et al. 9] , all of which claim 
a triplet pairing, are also explained in terms of the singlet 
scenario by Zutic and Mazin [1(3] and Mineev [ll|]. The 
most direct and virtually only probe to detect its parity 
is the KS experiment. In fact KS experiments using vari- 
ous nucleus, such as 87 Sr, 101 Ru, "Ru and 17 atoms, fail 
to pin down the spin direction of pairs, i.e. orientation 
of the (i-vector because of the invariance of KS for both 
field directions of c- and afr-axes as low as H — 200G [13] • 
There is no field direction where KS changes below T c . 
Thus at present it is fair to say that the two scenarios 
either based on singlet and triplet pairings are still under 
debate. Note that the appearance of magnetic field below 
T c associated with spontaneous time reversal symmetry 
breaking observed by /iSR experiment 13[ is explained 



equally by spin singlet scenario as well as triplet one 14 1 . 

We examine the parity issue in Sr2Ru04 through anal- 
yses of the specific heat experiment by Deguchi et al. [l6| 
under various T and H. There are several outstanding- 
problems posed by this experiment, whose understand- 
ing leads to a new clue for this debate. One of the most 
interesting discoveries is why the field dependence of the 
Sommerfeld coefficient j(H) = HmT^oC/T (C is the spe- 
cific heat) in the basal plane shows a concave curvature 
in spite of the existence of the line node gap. Namely, 
this is quite at odd because j(H) is expected to be a 
like behavior with a convex curvature due to line nodes, 
i.e. the so-called Volovik effect flit ]. It is remarkable to 
see that the concave curve becomes a Volovik s/H curve 
with a convex curvature when the direction of the ap- 
plied field moves away only by a few degrees of angle 9 
from the basal a6-plane (see inset (a) in Fi g. 3) . In ad- 
dition to analyses of the specific heat data[l6( we also 
examine magnetization datafl7j at low temperatures un- 
der a field. We explain these experiments based on an 
idea that strong Pauli paramagnetic effect is important 
in the basal ab plane physics of Sr 2 Ru04 and establish a 
consistent picture for its superconductivity. 

We calculate the vortex lattice state properties by qua- 
siclassical Eilenberger theory in the clean limit [HI ■ This 
framework is valid when kp£ 3> 1 (kp Fermi wave num- 
ber and £ coherent length), which is satisfied by Sr2Ru04. 
We include the paramagnetic effects due to the Zee- 
man term ^B-B(r). The flux density of the internal field 
is B(r) and /ie is a renormalized Bohr magneton [l9j |. 
The quasiclassical Green's functions g(uji + i/i.B,k, r), 
f(u>i + ifiB, k, r) and P(ui + ifiB, k, r) are calculated in 
the vortex lattice state by the Eilenberger equation 



{u n + ip,B + v(k F ) • [V 
{oj n + i/LB - v(k F ) • [V 



iA(r)]}/ = A(r)0, 
iA(r)]}/t = A*(r) ff , 



where g = (1 — ff^) 1 ^ 2 , Keg > 0, and the normal- 
ized Fermi velocity v is introduced so that (v 2 )k = 1 
where (• • -)k indicates the Fermi surface average. The 
paramagnetic parameter is jl = ^BQ/irk^T c . We con- 
sider the d-wave pairing for a pairing function with line 
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nodes on the two-dimensional (2D) cylindrical Fermi sur- 
face. The pair potential is selfconsistently calculated. 
The vector potential A for the internal magnetic field is 
selfconsistently determined by considering both the dia- 
magnetic contribution of supercurrent and the contribu- 
tion of the paramagnetic moment. We consider the large 
Ginzburg-Landau parameter k = 20. The local density 
of states is given by N(r,E) = N +1 (r,E) + N^x(r,E) 
with N a (r,E) = (Rc{g(uji + iafbB, k, r)\ iLU[ ^ E +i V })k for 
each spin component a = ±1. We typically use 77 = 0.01. 
The density of states is obtained by its spatial average as 
N(E) = {N(r, E)) r , which is identified as the Sommer- 
feld coefficient j(H) in specific heat at lower T. Using 
the Doria-Gubernatis-Rainer scaling, we calculate mag- 
netization M including diamagnetic and paramagnetic 
contributions. The details are found in Refs. (IR 

The paramagnetic parameter jl oc H°2 h /H p , which is 
a key parameter to analyze j{H), is related to the ra- 
tio of the hypothetical orbitally limited upper critical 
field H°l h and the Pauli limiting field H p = A /a/2/x b 
(Ao is the gap amplitude at T = 0). H p is a material- 
specific bulk parameter independent of the field orien- 
tation evidenced by nearly isotropic bulk susceptibility 
observed @ . The angle-dependence of the paramagnetic 
parameter ji(9) comes through the factor: H°^°{9). This 
orbital- limited H°^°{9) is sensitive to the field orienta- 
tion for highly anisotropic system such as in the present 
layered material; Sr 2 Ru04. 

The reduction of H c2 from H°^° due to the paramag- 
netic effect is obtained by solving the Eilenberger equa- 
tion as H c2 {p) = ff°2 rb /\A + 2.4/i 2 . This is derived orig- 
inally in dirty limit s-wave case plj . but we confirm it 
to be valid numerically in the present clean limit d-wave 
case too as seen from Fig. 1 where the calculated values 
are compared with this expression. 

It is natural to consider that H°^ a (9) is described 
by the effective mass model, namely H°2 h (9) / H° 2 ^ ab = 

l/\/r 2 sin 2 9 + cos 2 9 which simply embodies the fact 
that the orbital motion of electrons is determined by the 
directional cosine of the field to the basal plane. The 
anisotropy T = H°^ ab / H°^ c is an unknown parameter 
here. But it is assigned by the requirement that the ex- 
perimental H c2 (9) be reproduced theoretically. Namely, 
once T is determined, the angle dependence of H c2 (9) 
is automatically known through the angle dependence 
of the paramagnetic parameter (1(0), which controls the 
reduction of the upper critical field H c2 from the "hypo- 
thetical" orbital-limited field ff° 2 rb . 

Having known the paramagnetic depairing effect on 
H c2 iff) , we can calculate the angle dependence of the ob- 
served H c2 (9) where we take account of the fact that jl oc 
H c2°/ H p is ^-dependent through the factor H°^(9) given 
above. Thus we obtain ji(9) = jlo/\/T 2 sin 2 + cos 2 9 
with jlo being the value at 6 = 0. By combining these 
relations, we finally obtain the 9 dependence of the ob- 




FIG. 1: Reduction of H C 2(jx) as a function of fi evaluated 
by quasiclassical Eilenberger equation. The fitting curve is 
described well by H c2 (jl)/H°^ = 1/^1 + 2.4/P. 



served H c2 {9) as H c2 {9) = l/\/V 2 sin 2 9 + cos 2 9 + 2Ap, . 
This takes account of both orbital- and paramagnetic de- 
pairing effects simultaneously. In order to reproduce the 
observed anisotropy T ohs = 20, we find jlo = 3.41 when 
r = 107. Note that jlo and T are not independent pa- 
rameters. As shown in Fig. 2 our effective mass model 
with the paramagnetic effect explains the angle depen- 
dence of H c2 (9) once we fix one adjustable parameter. It 
is to be noted as shown in inset of Fig. 2 the jl(9) value 
is completely determined by the effective mass form with 
r = 107. 

As for the assigned V = 107 we point out that the dia- 
magnetic orbital current is determined by the perpendic- 
ular component of the average Fermi velocity to the field 
direction. Thus T is the anisotropy ratio of the Fermi 
velocities, namely T = \J (vp^ c ) / (vp ^ ab ). This quantity 

is determined directly by dHvA experiment; T a = 117, 
= 57 and T 7 = 174 for three bands a, (3 and 7 
respectively 5]. Note that a simple geometric average 
r o ff = ^(T a + + T^) = 116 is well compared with 
our assignment T — 107. In this sense there is virtually 
no adjustable parameter in our analysis. In passing we 
note that the observed ratio r obs = H c2 \\ ab / H c2 \\ c = 20 is 
strongly reduced from r o ff, apparently suggesting some 
reduction mechanism. We clarified it here. 

Let us now come to our main discussions on the anal- 
yses of the specific heat at a low T. In Fig. 3 we display 
7 (J?) for several values of jl together with the experi- 
mental data in inset (a) for various 9 values. They show 
strikingly similar behaviors as a whole. The larger angle 
data exhibit a strong upward curvature, corresponding 
to the conventional j(H) ~ \J~H which is characteristic 
to the line node gap structure. Those are reproduced 
in our /i=0.02, or 0.41 curves. As 9 becomes smaller, 
this changes into almost linear or concaved curves near 
H c2 . This behavior is captured by the theoretical calcu- 
lations for larger jtt's. Thus the overall "metamorphosis" 
of j(H) from the conventional yH to a strong convex 
curve is reproduced by increasing jl. As shown in inset 
(b) of Fig. 3, the data are fitted well by our calculations 
near H c2 where we have used the jl{9) values determined 
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FIG. 2: (color online) Calculated angle dependence of H C 2{9) 
(solid line). Circles [161 ] (squares [201]) are experimental data. 
Enlarged figure is shown in inset for small angles. The dotted 
line is the original orbital limit H°^°(6) of the effective mass 
form with V = 107. The dotted line also shows fi(6) with 
fo = 3.41 (right hand scale). 9 is the angle from the ab plane. 

above (see the inset of Fig. 2 with // n = 3.41). We have 
computed the six cases shown in Fig. 3 for fi values and 
obtained ^(H) for other /Ts by interpolation. 




0.4 H/Hc2 0.6 

FIG. 3: (color online) Zero-energy DOS j(H) at T = 0.1T C 
for ft = 0.02, 0.41, 0.86, 1.71, 2.57 and 3.41 from top to 
bottom. Inset (a) shows the experimental data [l(| for #=0° , 
2.5°, 3.0°, 5.0° and 90° from bottom to top. Inset (b) is 
the fitting of the data 0=0° by £=3.41, 0.5° (/i=2.36), 5° 
(/i=0.33) and 90° (/i=0.03) from bottom to top, which are 
shifted upwards. 




FIG. 4: (color online) (a) -y(H) for ft = 3.41, 0.60, 0.36,0.18 
and 0.06 from bottom to top. (b) Corresponding data [16J for 
6»=0°, 3°, 5°, 10° and 30°. 



In Fig. 4 we display the theoretical j(H) behaviors 
(a) and the corresponding specific heat data [lj| (b), 
where we read off fi(9) from the inset of Fig. 2. Our the- 
oretical curves explain these data in a consistent man- 
ner. In particular, it is noteworthy; (1) At 9=0° where 
fi(0) = /io = 3.41 is largest, ^(H) shows a y/H-\ike sharp 
rise in smaller H region because of the presence of line 
nodes. But it is limited only to lower fields. (2) In the 
intermediate wide field region (0.5T< H <1T), *y(H) ex- 
hibits an almost linear change in H. This extended linear 
change is shown to be consistent thermodynamically with 
magnetization M(T,H) behavior as seen shortly. (3) In 
the high field region (H >1T) towards H c2 = 1.5T, j(H) 
displays a sharp rise with a strong concave curvature. As 
H increases, the Pauli effect proportional linearly to H 
becomes growingly effective, modifying 'y(H) from usual 
\f~H to a concave H a -\\ke curve with a > 1. 

The data for 9=3° where fi{9 = 3°) = 0.60 show a sim- 
ilar behavior to that at 9=0°, but the features associated 
with the Pauli effect, namely, the existence of the inflec- 
tion point from convex to concave curves and sharp rise 
towards H C 2 are weaken. The j(H) data for higher an- 
gles (9 > 3°) exhibit an intermediate behavior between 
those at 8 = 0° and the ordinary \f~H curve, continu- 
ously changing its shape with 8. It is remarkable that 
the strong concaved curves of the experimental data for 
small angles, which were unexplained before, are repro- 
duced by the Pauli paramagnetic effect. Physically, this 
effect makes the conventional Abrikosov vortex state un- 
stable, ultimately leading to the normal state via a first 
order transition or the FFLO state. The sharp rise in 
'y(H) near H C 2 is a precursor to it. 

In Fig. 5 we show the calculated results of magneti- 
zation M(H) for several T's (a) together with the ex- 
perimental data[17j (b) to qualitatively understand the 
paramagnetic effects on M(T, H). We do not attempt to 
reproduce the data quantitatively because the data are 
in a qualitative nature due to hysteresis effects. It is seen 
from Fig. 5(a) that the magnetization with a convex cur- 
vature at lower field changes into that with a concave one 
towards H C 2 . There is an inflection point field Hk in be- 
tween. The relative position of Hk to each H C 2 decreases 
with T (also see insets). In higher T's Hk becomes in- 
visible because of thermal effect. These two features are 
observed experimentally as seen from Fig. 5(b). The 
inflection point field Hk roughly coincides with that in 
'j(H) as seen from Fig. 4, implying that these are ther- 
modynamically related to each other. 

As is seen from Fig. 5 upon lowering T the 
slope of M{H) at H c2 becomes steeper, meaning 
that K2 decreases, instead of increases as in usual 
superconductors [2lj. This is another obvious support- 
ing evidence that the paramagnetic effect is important in 

Sl'2Ru04. 

It is easy to derive a thermodynamic Maxwell rela- 
tion -jfj^ = -gjvM(T,H) from which we can see at low 
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FIG. 5: (color online) (a) Calculated magnetization curves 
for various T/T c = 0.1, 0.2, 0.3 • • • , 0.9 from bottom to top for 
jl — 1.71. Inset shows H C 2 and the inflection point Hk- (b) 
Corresponding data 17] for T/T c =0.1, 0.28, 0.40 and 0.56 
from bottom to top for H || ab. Inset shows H C 2 and "kink" 
field Hk in their terminology [17J ■ Magnetization of the nor- 
mal paramagnetic moment is substracted. 



T, = /3(H) with M(T,H) = M (H) + \P{HW 2 - 

We estimate (3(H) from the experimental data [Ljj in 
Fig. 5, finding that (5(H) ~ const for 0.5T<£T<1T 
and 0(H) cx H 3 for lT<i/<1.35T. This implies that 
-y(H) cx H(H*) for 0.5T<H<1T (lT<i?<1.35T). These 
behaviors in j(H) are indeed seen for the 6 = 0° data 
shown in Fig. 4. These analyses, which are free from any 
microscopic model, mean that the mysterious behavior 
of j(H) is supported to be true thermodynamically and 
comes from the intrinsic nature deeply rooted to the su- 
perconductivity in Sr2Ru04. 

There are several known difficulties associated with 
the most popular two component chiral p-wave pairing; 
z(p x + ip y )[22] or z(p x + ip y ) cosp z [23j: Experimentally 
these triplet states are unable to explain the paramag- 
netic effects mentioned above because the d- vector is not 
locked in the basal plane. Theoretically these states give 
a large in-plain H C 2 anisotropy 24j which is not observed. 
The present singlet scenario is free from it. 

Let us go on considering the high field phase for H\\ab 
observed as the double transition[20(. It appears in a 
narrow H-T region along H c2 \\ a b, starting at T = 0.8K, 
or To = 0.53T C at which three transition lines meet, 
giving rise to a tricritical point in H vs.. T plane. 
To is remarkably similar to the so-called Lifshitz point 
Tl = 0.56T C in the FFLO phase diagram for a Pauli 
limited superconductor where the orbital depairing is 
quenched completely. This number Tl = 0.56T C is uni- 
versal, valid for a variety of situations, including 3D 
Fermi sphere s- wave [25$. 2D s-wave[26j and <j-wave|27j. 
and ID s-wavepjj models. Our identified large param- 
agnetic parameter jl = 3.41 means that our system is 
in almost Pauli limiting where the orbital effect is al- 
most perfectly quenched because the two-dimensionality 
in Sr2Ru04 is so extreme. In fact note that the iden- 
tified anisotropy r=107 implies H°$ ab ~ 7.5T which is 



reduced to H c2 \\ a i, = 1.5T by the Pauli effect. Thus we 
propose here to identify this high field phase as FFLO. 

The extreme two-dimensionality is obvious: If H is 
tilted away from the ab plane only by 9 > 0.3°, the 
double transition vanishes (20j. According to Nakai, et 
al.^ the FFLO region at low T occupies - 0.8% be- 
low H C 2, which is comparable with the width ~200G 
of the high field phase below H C 2\\ a b = L5T, a re- 
gion 200G/1.5T~1.3%[2^. Guided by the known phase 
diagram]^, we predict that as the field orientation 9 in- 
creases, jl decreasing, this high field phase survives only 
for < 8 < 0.3° and quickly diminishes for 9 > 0.3°. 
At around 9 ~ 1.0° there appears a first order transition 
along H c2 line instead of FFLO. Then for 9 > 2.0° it also 
disappears above which the paramagnetic effect becomes 
ineffective and Sr2Ru04 is described by a conventional 
singlet superconductor with line nodes. These predic- 
tions based on our analyses are all testable experimen- 
tally although the details should be further sharpened 
theoretically. 

In conclusion, we have analyzed both specific heat at 
lower T and magnetization M(T, H) by self-consistently 
solving microscopic quasi-classical Eilenberger equation 
for the gap function with line nodes. It is seen that the 
Pauli paramagnetic depairing effect is essential in under- 
standing the data in Sr2Ru02- This is possible only for 
either singlet pairing, or triplet pairing with the d vector 
locked in the basal plane. 

We are grateful for useful discussions and communi- 
cations to K. Deguchi, K. Tenya, K. Ishida, Y. Maeno, 
H. Adachi, N. Nakai, P. Miranovic and Y. Matsuda. M. 
Tsutsumi helps us for preparing figures. 
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